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The Swift-BAT hard X-ray sky monitoring unveils the orbital period 
of the HMXB IGR J16493-4348 
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ABSTRACT 

IGR J 16493^1348 is a supergiant high mass X-ray binary discovered by INTEGRAL in 2004. 
The source is detected at a significance level of ~ 21 standard deviations in the Swift-BAT 
survey data collected during the first 54 months of the Swift mission. The timing analysis 
reveals an orbital period of ^6.78 days and the presence of a full eclipse of the compact object. 
The dynamical range (variability up to a factor ^20) observed during the BAT monitoring 
suggests that IGR J 16493^-348 is a wind-fed system. The derived semi-major axis of the 
binary system is ~ 55i?Q with an orbit eccentricity lower than 0.15. 
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1 INTRODUCTION 



The IBIS/ISGRI telescope Jubertini et all 120031: iLebrun et all 
12003b on board the INTEGRAL satellite JWinkler et alj|2003l) has 
detected a large number of new supergiant high mass X-ray binaries 
(sgHMXB) characterized by a high absorption column (Nh > 10 23 
cm 2 ) and/or by fast bright transient events. This result was achieved 
thanks to a continuous monitoring of the Galactic plane with deep 
exposure of the Galactic centre and thanks to a combination of IS- 
GRI large field of view, good sensitivity, and wide energy range. 
The fourth IBIS/ISGRI survey catalogue ( Birdet alj|20ld) reports 
the discovery of 32 new objects that have been classified as HMXB 
and 84 unidentified sources with a Galactic latitude \b\ < 5°. 

The Bu rst Alert Telescope (BAT. iBarthelmv et alj|2005h on 
board Swift (Gehrels et al. 2004) has been performing a continu- 
ous monitoring of the sky in the hard X-ray energy range (15-150 
keV) since November 2004. The telescope, thanks to its large field 
of view (1.4 steradian half coded) and its pointing strategy, covers 
a fraction between 50% and 80% of the sky every day. This has 
allow ed the detection of m any of the new INTEGRAL HMXBs 
(e.g. ICusumano et al.ll2010r > and the collection of their long term 
light curves and of their spectral energy distributions. The long and 
continuous monitoring of these sources allows to investigate the in- 
trinsic emission variability, to search for long periodicities (orbital 
periods) and to discover the presence of eclipse events. The role 
of Swift-BAT is therefore fundamental to unveil the nature and the 
geometry of these binary systems. 

In this Letter we analyze the soft and hard X-ray data col- 
lected by Swift on IGR J 16493-4348. Th is source was discovered 
by INTEGRAL in 2004 JBird et all2004l) and it was initially asso- 
ciated with the radio pulsar PSR J1649-4349 because of a spatial 



coincidence. A later INTEGRAL observation with a deep expo- 
sure allowed to re duce the positional un certainty and to reject the 
pulsar association JGrebenev et alj|2005l) . A follow-up observation 
with Chandra found a soft X-ray counterp art at RA(I2000) = 16 fe 
49 m 26.92 s ; Dec(I2000) = -43° 49' 8.96" JKuiper et alj|2005l) al- 
lowing the optical a ssociation with 2MASS 11642695-43490 90, a 
B0. 5 lb supergiant JNespoli. Fabregat. & Mennickentl 120081) . The 
Chandra observation also revealed strong evidence for variabil- 
ity. The IGR 116493 -4348 energy distribution extracte d from two 
RXTE observations JMarkwardt, Swank.& S mith 2005) was mod- 
eled with a highly absorbed power law (Nh ~ 10 23 cm -2 ) with 
a photon index of ~1.4 and a 20-40 keV flux of 2.1 x 10~ n erg 
cm -2 s _1 . A spectral analysis using non simultaneous data from 
Swift-XRT and INTEGRAL (Hill et al. 2008) showed the presence 
of a cutoff at ~ 15 keV. The information derived from the X-ray 
observations and the identification of the spectral type of the opti- 
cal counterpart allowed the classification of IGR J16493-4348 as a 
supergiant High Mass X-ray Binary (sgHMXB). 

This Letter is organized as follows. Section 2 describes the 
data reduction. Section 3 reports on the timing analysis and in Sect. 
4 we discuss our results. 



2 DATA REDUCTION 

The raw BAT survey data of the first 54 months of the Swift mis- 
sion were retrieved from the HEA SARC public archiv qj and pro- 
cessed with a dedicated software JSegreto et al.l 120101) , that per- 
forms screening, mosaicking and source detection on BAT data 



http://heasarc.gsfc.nasa.gov/docs/archive.html 
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Figure 1. Top panel: 15-50 keV significance map in the neighborhood of 
IGR J16493— 4348. The color-bar represents the significance levels. 



and produces spectra and light curves for any given sky position. 
Figure Q] shows the 15-50 keV significance sky map (exposure 
time of 17.7 Ms, duty cycle ~ 13%) centered in the direction of 
IGR J16493— 4348. The source is clearly detected at a significance 
of 20.8 standard deviations. The light curve of IGR 116493-^348 
was extracted in the 15-50 keV energy range with the maximum 
available time resolution (~ 300 s) , subtracting the con tamination 
of the nearby sources as detailed in lSegreto et al.l ( 120100 . The time 
tag of each bin, defined as the central time of the bin interval, was 
corrected to the Solar system barycentre (SSB) by using the task 
EARTH2SUN. 

W e also re-analyzed the data of the Swift-XKT JBurrows et al.l 
120051) IGR J 16493-4348 observation performed on 2006 March 1 1 
(Obsld 00030379002), for a total exposure time of 5.6 ks. The data 
were processed with standard procedures (XRTPIPELINE v.0.12.4), 
filtering and screening criteria, using FTOOLS in the HEASOFT 
packa ge (v 6.8). The source was observed in Photon Counting 
mode JHill e t al. 2004) that provides a time resolution of 2.5 s. 

We adopted standard grade filtering 0-12. The source events 
were extracted from a circular region of 20 pixels radius (1 
pixel=2.36") centered on the source position as determined with 
XRTCENTROID. The event arrival times were converted to the SSB 
with the task BARYCORR. 



3 TIMING ANALYSIS 

3.1 BAT survey data 

The average count rate in the BAT light curve is 1.03 x 10 -4 count 
s _1 . When considering the light curve at the highest resolution the 
maximum of the deviation from the average rate is about 7a, cor- 
responding to an increase in rate of a factor ~ 20. 

We analyzed the long term BAT light curve to search 
for intensity modulations by applying a folding technique and 
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Figure 2. a: Periodogram of Swift-BAT (15-50 keV) data for IGR J16493- 
4348. Po, Pi, P2 and P3 are defined in Section 3. The inset shows a close- 
up view of the x 2 distribution around Po. b: Light curve folded at a period 
Po = 6.783 days, with 16 phase bins, c: Phase coverage of the light curve 
bins. The errors are the square root of the number of bins in each phase 
interval d: Distribution of the z=x 2 -F x (P) values extracted in the period 
range between 0.5 and 10 days excluding the z values obtained at Po, Pi, 
P2 and P3 . The continuous line is the best fit obtained with an exponential 
model applied to the tail of the distribution (z>25). 
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searching in the 0.5-1 00 d period range. The period resolution 
dBuccheri & Saccoll 19851) is given by P 2 /(N AT), where N = 16 
is the number of trial profile phase bins, and AT =140,213,639.0 s 
is the data time span. The average rate in each phase bin was evalu- 
ated by weighting the light curve rates by the inverse square of the 
corresponding statistical error 
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where Rj is the average rate in the j-th phase bin of the trial profile, 
Ti are the rate of the light curve whose phase fall into the j-th phase 
bin and eri are the corresponding statistical errors. The error on Rj 

is (yX/ V er f)~ ■ The weighting procedure was adopted to deal 
with the large span of er; and it is justified because the BAT data 
are background dominated. 

Figure [2] (a) shows the periodogram where several features 
emerge. The highest feature, with a \ value of ~ 240, is Po = 
6.782 ± 0.005 d, where the error is evaluated as the standard devi- 
ation of the gaussian that best fits the x 2 values around the Po peak 
with x 2 > 150. We also see other evident features at periods multi- 
ple of Po (2P , 3Po, 4P , 5P in Fig.H}. The pulsed profile (Fig.|2] 
b) folded at P with T opoch =54173.757 MJD, shows a flat intensity 
level, abruptly broken by a deep full eclipse. The phase coverage of 
the light curve bins (Fig. [2] c) shows that this dip is not due to an 
accidental under-sampling of the light curve at these phases. The 
centroid of the eclipse, evaluated by fitting the data around the dip 
with a Gaussian model, is at phase 0.319 ± 0.015 corresponding to 
Todipso = (54175.92 ± 0.10) ± nxP„ MJD. 

As a consequence of the time variability of the source and of 
the presence of a periodic signal, the average x 2 m the periodogram 
is far from the average value expected for white noise (N — 1) and 
the x 2 statistics cannot be applied to evaluate the significance of the 
detected periodicity. Therefore, we applied two alternative methods 
to evaluate the significance of Po. 

(i) We fit the periodogram with a second order polynomial 
[F X (P)] and subtracted the trend from the x 2 distribution. The value 
of z=x 2 -F x (P) at Po is 196.6. We therefore built the histogram of 
the z distribution (Fig. [2] d) from 0.5 to 10 d excluding the inter- 
val around Po and those around Pi, P2, and P3 as marked in Fig. [2] 
(a). The latter features are not noise fluctuations, but are tied to Po: 
Pj (=Po/2) is due to the presence of the deep eclipse, still visi- 
ble in the light curve folded at half of the period; P2 (~ 2.27d= 
(1/Po + I/Pi)" 1 ) and P 3 (~ 4.54d= [(1/P + l/Pi)/^ 1 ) are 
due to beat frequencies between Po and Pi. We fit the resulting dis- 
tribution for z> 25 with an exponential function and evaluated the 
integral of the best-fit function beyond z=196.6. This integral yields 
a number of chance occurrences due to noise of 4.6 x 10~ 1S , cor- 
responding to a significance for the detected feature of ~ 8.5 stan- 
dard deviations in Gaussian statistics. 

(ii) We generated 1000 light curves scrambling the observed 
rates while keeeping unchanged the temporal distribution of the 
bins. For each of them we have produced the periodogram in the 
0.5-10 d time interval (49335 trial periods). The highest value ob- 
tained among all these periodograms is 97.9. Therefore, the proba- 
bility of random occurrence for the observed x 2 value (~ 240) is 
lower than 2 x 10 -8 that corresponds to a significance of Po higher 
than 5.5 standard deviations. 

Figure [3] shows the 15-50 keV light curve of 
IGR J 16493-4348, obtained excluding all pointings where 
the source was at an off-axis angle higher than 40 deg, with a bin 
time of Po = 6.78 d, excluding bins with an exposure fraction less 
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Figure 3. 15-50 keV BAT light curve of IGRJ16493-4348 with a bin time 
Pq = 6.78 d. The dashed line represents the average count rate. 



than 5%. The average intensity in each Po bin shows a maximum 
variability of a factor ~ 4. 



3.2 XRT data 

The knowledge of the orbital period and of the middle eclipse time 
allows us to determine that the Swift-XKT observation was per- 
formed at the orbital phase interval 0.45-0.51 with respect to the 
epoch of the eclipse. The 0.2-10 keV light curve shows a persistent 
emission with an average count rate of 0.376±0.015 count s _1 and 
a variability within a factor 3. 

We applied the folding technique to search for short periodic- 
ities in the time range 5-1000 s. The XRT SSB arrival times were 
corrected for the pulsar binary motion using Po, T ec ii pse , a semi- 
major axis a ~ 55 Rq ~ 255, 365 lt-s (see Sect. 4 for its deriva- 
tion), and assuming an orbit inclination of 90 degrees and an ac- 
centricity e=0. We found no evidence for any significant periodic 
signal. 



4 DISCUSSION 

We have presented the results obtained from the analysis of the 
data collected by Swift-BAT during the first 54 months of the Swift 
mission on the supergiant HMXB IGR J 16493-4348. The source is 
detected at a significance level of ~21 standard deviations. 

The long monitoring coupled with a good observation duty 
cycle (~ 13% per day) allowed us to unveil a periodicity of ~ 
6.782 ± 0.002 d that we interpreted as the orbital period of the 
binary system. 

The knowledge of the orbital period allows us to derive the 
the semi-major axis of the binary system through the Kepler's third 
law: a 3 = G P(j(M* + M x )/4tt 2 where M* and M x are the 
masses of the supe rgiant and compac t object, respectively. Adopt- 
ing M* = 47 M & JSearle et alj2008l , for an B0.5 lb star, of radius 
Ri, — 32.2 Rq) and a typical neutron star mass of M% = 1.4 Mq 
we obtain a ~ 55 Rq ~ 2 R *. This orbital separation is com- 
mon among classical sgHMXB dWalter & Zurita Herasll2007l) . The 
folded light curve shows the presence of a full eclipse with an esti- 
mated time of mid-eclispe at T cc h psc = (54175.92 ± 0.10) ± nx 
Po MJD and it appears flat outside the eclipse. The duration of the 
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Figure 4. The Lagrangian point of the binary system as a function of the 
orbital phase for different orbital eccentricity. The straight dashed line rep- 
resents the radius of the companion star. 
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author. 



eclipse (~ 12% of the orbital period) is roughly consistent with the 
occultation of the compact source by the companion star in a highly 
inclined orbit with low eccentricity. 

The 15-50 keV light curve of IGR J 16493^1348 is variable 
up to a factor ~20 on timescales of ~300 s while its average over 
an orbital period shows a variability up to a factor ~4. The inferred 
orbital separation suggests a wind-fed system. This hypothesis and 
the radius of the supergiant allows to limit the eccentricity of the 
orbit. Figure [4] shows how the Lagrangian point between the two 
stars changes as a function of t he orbital phase for different orbit 
eccentricities ( IPacz vriskil 1 1 97 11) : the maximum eccentricity com- 
patible with a wind-fed system (Lagrangian point lower than the 
supergiant radius) is 0.15. 
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